Redshift space correlations and scale-dependent stochastic biasing of density peaks 
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We calculate the redshift space correlation function and the power spectrum of density peaks of a 
Gaussian random field. Our derivation, which is valid on linear scales k < 0.1 /iMpc -1 , is based on 
the peak biasing relation given in Desjacques [Phys. Rev. D. , 78, 3503 (2008)]. In linear theory, 
the redshift space power spectrum is 

Ppk(k, fJ.) = exp(-fa^ el k 2 ^ 2 ) [ftpk(fc) + bvei(fc) /M 2 ] 2 Ps(k), 

where fi is the angle with respect to the line of sight, cr vc i is the one-dimensional velocity dispersion, 
/ is the growth rate, and b p k(fc) and b VB i(k) are fc-dependent linear spatial and velocity bias factors. 
For peaks, the value of cr vc i depends upon the functional form of b vc \. When the fc-dependence 
is absent from the square brackets and 6 vc i is set to unity, the resulting expression is assumed to 
describe models where the bias is linear and deterministic, but the velocities are unbiased. The 
peaks model is remarkable because it has unbiased velocities in this same sense - peak motions are 
driven by dark matter flows - but, in order to achieve this, b vc i is k— dependent. We speculate that 
this is true in general: fc-dependence of the spatial bias will lead to fc-dependence of fo V ei even if the 
biased tracers flow with the dark matter. Because of the fc-dependence of the linear bias parameters, 
standard manipulations applied to the peak model will lead to fc-dependent estimates of the growth 
factor that could erroneously be interpreted as a signature of modified dark energy or gravity. We 
use the Fisher formalism to show that the constraint on the growth rate / is degraded by a factor of 
two if one allows for a fc-dependent velocity bias of the peak type. Our analysis also demonstrates 
that the Gaussian smoothing term is part and parcel of linear theory. We discuss a simple estimate 
of nonlinear evolution and illustrate the effect of the peak bias on the redshift space multipoles. For 
k < 0.1 /iMpc -1 , the peak bias is deterministic but fc-dependent, so the configuration space bias 
is stochastic and scale dependent, both in real and redshift space. We provide expressions for this 
stochasticity and its evolution. 
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I. INTRODUCTION 



While velocities are directly measured through their 
Doppler (red)shifts, accurate measurement of cosmologi- 
cal distances are only available for nearby cosmic objects, 
and even at these small scales they are plagued with ob- 
servational biases. Therefore, most observational data is 
described in terms of redshifts, e.g. three-dimensional 
(3D) galaxy surveys provide the angular positions and 
redshifts of galaxies. Redshifts differ from distances by 
the peculiar velocities (deviations from pure Hubble flow) 
along the line of sight. These generate systematic differ- 
ences between the spatial distribution of data in redshift 
and distance (or real) space which are commonly referred 
to as redshift distortions Kaiser [2[ first derived an 
expression which describes the effect of linear peculiar 
motions on 3D power spectra. References Q and [H pro- 
vide two very different derivations of this same expres- 
sion. Whereas the original derivation made no assump- 
tion about the form of the density and velocity fields, the 
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other two assume they are Gaussian distributed. 

The Kaiser formula has been used to interpret observa- 
tions of the redshift space clustering of galaxies. The an- 
gular dependence of the redshift distortions can be used 
to measure the logarithmic derivative / = dlnZ?/dlna or 
growth rate @ at multiple redshifts and thus potentially 
constrain many of the dark energy or modified gravity 
models (e.g. Q; for a review of these scenarios, see @). 
Essentially all analyses to date assume that i) galaxies 
are biased tracers of the underlying matter field, ii) the 
bias is linear, local and deterministic [2|, y, l8Hl0l| and iii) 
the velocities of the tracers are unbiased. In fact, ex- 
cept on the largest scales, the relation between the dark 
matter and galaxy fields is almost certainly nonlinear, 
nonlocal, and scale dependent [Tl| . Our main goal in the 
present study is to explore what complexities one might 
expect on smaller scales where the bias relation is more 
complicated, and where the velocities may also be biased. 
We do so by investigating the impact of redshift distor- 
tions on the correlation function of density maxima in a 
Gaussian density field. 

We have chosen to study density peaks because the 
statistics of Gaussian random density 12[ and velocity 
fields [l3] in a cosmological context, and of the peak dis- 
tribution in particular, has already received considerable 
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attention [T3 - l20j . Some of these results have been used 
in studies of the nonspherical formation of large-scale 
structures piTt^ . Others, especially from peaks theory, 
have been used to interpret the abundance and cluster- 
ing of rich clusters [25l428l | . Density peaks define a well- 
behaved point-process which can account for the discrete 
nature of dark matter halos and galaxies. On asymp- 
totically large scales, peaks are linearly biased tracers of 
the dark matter field, and this bias is scale independent 
[3 E3, [HI- However, these conclusions are based on a 
configuration-space argument known as the peak back- 
ground split. Extending the description of peak bias to 
smaller scales is more easily accomplished by working in 
Fourier space. It has been shown that peaks are lin- 
early biased with respect to the mass, but this bias is 
/c-dependent pH.[29j. 

The first part of this paper demonstrates that, in the 
large-scale limit, the configuration and Fourier-based ap- 
proaches yield consistent results. This is important, be- 
cause the first (and only other) study of the redshift space 
clustering of peaks, reference [20T ]. reported that in red- 
shift space peaks behave very differently from the de- 
terministic, linear and scale independent biased tracers 
investigated in 0, [H H, [l(| ■ Since the linear bias as- 
sumption that is extensively advocated to convert large 
scale redshift space measurements into information about 
the background cosmology d, [3l|, the fact that peaks 
might behave very differently is potentially very worry- 
ing. In addition, peak velocities exhibit a fc-dependent 
bias even though peaks locally flow with the dark matter 
29] . This is remarkable given that one commonly refers 
to such flows as having unbiased velocities. We explain 
the origin of this effect and argue that it should be a 
generic feature of any fc-dependent spatial bias model. 
Again, however, peaks are remarkable because, in the 
high peak limit where their spatial bias is expected to be 
linear and scale independent, their velocity bias remains 
fc-dependent. 

The second part shows that, at the linear order, red- 
shift space distortions for peaks can be recast in a way 
that retains the simplicity of the original Kaiser formulae 
while generalizing them to tracers whose linear bias 
is fc-dependent. Because the present derivation is based 
on a model which is supposed to be accurate at smaller 
scales, we can identify an important term which does not 
appear in Furthermore, our analysis reaches very 
different conclusions from that of [20| . Our peaks-based 
formula for redshift space distortions, which includes fc- 
dependent linear bias factors for both the density and 
the velocity fields, has a rich structure. We hope it will 
serve as a guide for what one might expect in the case of 
more realistic (nonlinear, nonlocal, scale dependent) bias 
prescriptions. 

In the last part of this study, we use the Fisher for- 
malism to quantify the extent to which any fc-dependent 
velocity bias of the peak type would degrade the uncer- 
tainties on the growth rate /. We also demonstrate the 
stochastic nature of the peak bias and discuss its evolu- 



tion with redshift. The peak biasing is interesting be- 
cause, although it is deterministic in Fourier space, it is 
stochastic in real space. A final section summarizes our 
findings and speculate on some implications of the peak 
model. 

Throughout the paper we work in the "distant ob- 
server" limit, where the line of sight is oriented along 
the z direction. In all illustrative examples, we assume a 
flat ACDM cosmology with fl m = 0.279, fl h = 0.0462, 
h = 0.7, n s = 0.96 and a present-day normalisation 
as = 0.81 |34(. It will also be convenient to work with 
scaled velocities Vi = Vi/(aH /), where v^ is the (proper) 
peculiar velocity, H = dlna/dt, and / = d\nD/dh\.a with 
D(z) the linear theory growth factor. At z = 0.5 this is 
aHf ss 61 kms -1 Mpc -1 . As a result, m has dimen- 
sions of length. 

II. PROPERTIES OF DENSITY PEAKS 

We begin by reviewing some general properties of 
peaks in Gaussian random fields. We then discuss the 
biasing relation which is used in the calculation of the 
redshift space correlation of density maxima. 

A. Spectral moments 

The statistical properties of density peaks depend not 
only on the underlying density field, but also on its first 
and second derivatives. We are, therefore, interested in 
the linear (Gaussian) density field 5(x) and its first and 
second derivatives, diS(x) and didj5{x). In this regard, 
it is convenient to introduce the normalised variables 
v = 5{x)/<jq and u = —\7 2 5(x)/<j2, where the a n are 
the spectral moments of the matter power spectrum, 

1 r°° 

al ^— dk fc 2 (" +1 > P s (k, z)W(k, R S ) 2 . (1) 
2tt 2 Jo 

Here, Pg(k,z) denotes the dimensionless power spec- 
trum of the linear density field at redshift z, and W 
is a spherically symmetric smoothing kernel of length 
Rs (a Gaussian filter will be adopted throughout this 
paper) introduced to ensure convergence of all spectral 
moments. We will use the notation Ps s (fc, z) to denote 
Ps(k, z)W(k, Rs) 2 - The ratio a jo\ is proportional to 
the typical separation between zero-crossings of the den- 
sity field pj}. For subsequent use, we also define the 
spectral parameters 

In = (2) 

which reflect the range over which k 2n+1 P$ s (fc, z) is large. 

We will also need the analogous quantities to a 2 but 
for non-zero lag: 

i r°° 

^ (r)= 2^i dkk2{n+1)p Ss(k,z)3t(kr), (3) 
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where je(x) are spherical Bessel functions. As I gets 
larger, these harmonic transforms become increasingly 
sensitive to small-scale power. 

Finally, we note that the auto- and cross-correlations 
of the fields V{, 5, di8 and didjS can generally be de- 
composed into components with definite transformation 
properties under rotations. Reference [29j gives explicit 
expressions for the isotropic and homogeneous linear den- 
sity field. 



B. Smoothing scale and peak height 

The peak height v and the filtering radius Rs could 
in principle be treated as two independent variables. 
However, in order to make as much connection with 
dark matter halos (and, to a lesser extent, galaxies) as 
possible, we assume that density maxima with height 
v = S sc (z) / ao(Rs) identified in the primeval density field 
smoothed at scale Rs are related to dark matter ha- 
los of mass Ms collapsing at redshift z, where 5 sc (z) 
is the critical density for collapse at z in the spherical 
model [13, HH . For sake of illustration, we will present 
results at z = 0.5. In the background cosmology we as- 
sume, the linear critical density for (spherical) collapse 
at z = 0.5 is S sc w 1.681. The Gaussian smoothing scale 
at which v = 1 is R+ « 1.3 h~ 1 Mpc, so the characteristic 
mass scale is M* 6.5 x 10 11 M //i. 

While there is a direct correspondence between mas- 
sive halos in the evolved density field and the largest 
maxima of the initial density field, the extent to which 
galaxy-sized halos trace the initial density maxima is un- 
clear. Therefore, we will only consider mass scales Ms 
significantly larger than the characteristic mass for clus- 
tering, M*, for which the peak model is expected to work 
best. We will present results at redshift z = 0.5 for 
two (Gaussian) filtering lengths, Rs = 2.5 /i _1 Mpc and 
Rs = 4 /i _1 Mpc; these correspond to masses Ms = 1.9 x 
10 13 M Q /h and 7.8 x 10 13 M e /h, which roughly match 
the mean redshift and typical mass of halos harbouring 
luminous red galaxies (LRGs) in the Sloan Digital Sky 
Survey (SDSS) (H-H^l- This makes a /(Ti ~ 3.2 h^Mpc 
and 4.9 h~ 1 Mpc. To help set scales in the discus- 
sion which follows, the associated values of (i/,b^,b<^) 
are (2.1, 1.0, 16.4 /i 2 Mpc" 2 ) and (2.8, 2.8, 43.0 /i 2 Mpc~ 2 ). 
The three-dimensional velocity dispersion of these peaks 
is (T 2 _ 1 (1 — 7q): for our two smoothing scales, this cor- 
responds to (7.12 /i _1 Mpc) 2 and (6.66 /i _1 Mpc) 2 (recall 
that our velocities are in units of aH /, so they have di- 
mensions of (length) 2 ). 



C. Biasing 



Sn pk (x) = b v 8s{x) - b c V 2 S s (x) 
Vpk(x) = v s (x) |Vfe(x) , 



(4) 



where vs is the dark matter velocity smoothed at scale 
Rs (so as to retain only the large-scale, coherent motion 
of the peak) , and the bias parameters b v and b^ are 



-(- 

1 ( u — ~fll> 
02 



7i it 



CO 



(5) 



Here, u denotes the mean curvature of the peaks. Fur- 
thermore, b v is dimensionless, whereas b^ has units of 
(length) 2 . Note that b v is precisely the amplification fac- 
tor found by [l7| who neglected derivatives of the density 
correlation function (i.e. their analysis assumes 6^ = 0). 
We emphasize that Eq.Q is the only bias relation that 
can account for the first order peak correlation and mean 
streaming. 

Strictly speaking, the bias relation Q is nonlocal be- 
cause of the smoothing. In configuration space, the peak 
bias 6 p k at first order could thus be defined as the con- 
volution 



(b pk ® S) (x) = (b v - 6 C V 2 ) S s (x) , 



In Fourier space, this becomes 



b pk (k) = (b v + b c k 2 ) W(k, Rs) 



(6) 



(7) 



so it has the same functional form as Eq. (57) in reference 
[T]| who considered density extrema. Our coefficients 
thus agree with those of [Tl| only in the limit v 1, in 
which nearly all extrema are local maxima. 

This bias relation is distinct from either linear [l4j or 
nonlinear 38-40] biasing transformations of the density 
field for which 6^ — 0. Note in particular that Eq. (O 
shows that local bias schemes can generate fc-dependent 
bias factors if the bias relation involves differential oper- 
ators. Furthermore, when v ^> 1, then u — > 71^, so that 
co^i^ - > v and o^fr^ — > [30| . This is clearly seen in Fig.Q] 
where the biasing factors are plotted as a function of the 
peak height. Thus, the spatial bias of the highest peaks is 
expected to become scale independent, approaching the 
local deterministic relation of linearly biased tracers for 
which there is no fc-dependent bias. However, notice that 
the fc-dependence in the velocity bias remains. We will 
return to this point shortly. 



The large-scale asymptotics r — > 00 of the two-point 
correlation ^pkM and line of sight mean streaming [^12 • 
r] (r) for discrete local maxima of height v can be thought 
of as arising from the continuous, nonlinear bias relation 



D. Relation to peak background split 

There is another route for estimating large scale bias 
of peaks [26| which utilizes the peak background split 
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argument [17], I42H44I ] . This approach which is very dif- 
ferent from ours, because it is based on configuration 
space counts-in-cells statistics. In particular, it makes no 
mention of the bias in Fourier space. 

The large scale bias predicted by this approach is [26[ 



&pkbs = — 



1 d\n[n pk (v)] 



where 



'pk 



i 



d \nv 



^ 2/2 G (7i,7i^) 



(8) 



(9) 



(27r) 2 i?3 

is the differential averaged number density of peaks in 
the range v to v + dv [l3]- Here R\ — v / 3cri/o- 2 oc 
Rs characterises the typical radius of density maxima, 
and G is given by setting n = in our equation (jAlOl) . 
Therefore, 



frpkbs — 



v 2 + gi 

a v 



where 



91 



91nG (7i,2/) 



d\ny 

Performing the derivative yields 

Gi (71, 71 v) /G (71, 71^) - 71 v 



,9i 



-7!i/ ■ 



(10) 



(11) 



(12) 



where Gi is given by equation (|A10[) with n = 1. How- 
ever, G1/G0 = u (see the discussion immediately follow- 
ing equation I AlOp . so 



.91 = -71^ 



71 v 



7f 



= -71!/ (72 & 



(13) 



The last equality follows from the definition of (Eq. [5]) . 
Equations © and ([5]) eventually imply that 



71 f 02 &c = _ b v <Jo), 



so 



frpkbs — 



v 2 - v{y - bvOo) 



= K 



(T V 



(14) 



(15) 



This demonstrates that the large-scale, constant, deter- 
ministic bias factor returned by the peak background 
split approach is exactly the same as in our approach, 
when one considers scales large enough such that the In- 
dependence associated with the term can be ignored. 

This is very reassuring for two reasons. First, recall 
that our exp ressions for b v and b^ only agree with those 
given in [111 ] in the limit v 3> 1 (in which extrema are al- 
most certainly peaks). The analysis above shows that our 
b v is the appropriate generalization to lower v. And sec- 
ond, the peak background split approximation has been 
shown to provide an excellent description of large scale 
peak bias in simulations (2(| |4l| . Since our expressions 
reproduce this limit, we have confidence that our ap- 
proach will provide a good approximation on the smaller 
scales where the peak-background split fails (i.e., where 
the bias b pk becomes scale dependent). 



E. Power spectra and correlation functions 

Using the bias relations ((4]), it is straightforward to 
show that the real space cross- and auto-power spectrum 
are 



PpM(fc) = (b v + b c k 2 ) P s (k) W(k, R s ) 
Ppk(fc) = {K + b c k 2 ) 2 P Ss (k). 



(16) 
(17) 



We have omitted the explicit redshift and ^-dependence 
for brevity. The corresponding relations for the correla- 
tion functions are 



^Ar)-K^\r)+b c ^\r) 
U(r)=bltfXr)+2b„b c £\r) 



b 2 (r)$\r) 



(18) 



(19) 



where the final expression defines the (scale dependent) 
peak bias factor in configuration space. As shown in [29( 
(for £ p k(?")) and in Appendix fAl ffor £ p k, 5 (?")), these ex- 
pressions agree with those obtained from a rather lengthy 
derivation based on the peak constraint, which involves 
joint probability distributions of the density field and its 
derivatives. It is worth noticing that, while expressions 
(TT71) and (TT9"|) for the auto-power spectrum and correla- 
tion are only valid at first order in the correlation func- 
tions £p n ^ , the cross-power spectrum (|16[) and correlation 
(fT5|) are exact to all orders. 



F. Velocities 

In what follows, we will be interested in redshift space 
quantities, for which the velocity field also matters. The 
bias of peak velocities is particularly simple in Fourier 
space. Taking the divergence of Eq. (0}, we find 



9 pk (x) = V • v pk (x) = V • v s (x) - ^V z S s (x). (20) 

The linear continuity equation stipulates that 0s(x) = 
V • vs(x) = —5s(x), so the result of Fourier transforming 
the expression above implies that 



9 pk ( fc ) 



'0 7,2 



W(k,R s )6(k) = b vcl (k)0(k) 



(21) 

This defines the peak velocity bias factor, & ve i(fc), which 
depends on k but not on v. As seen in Fig. [TJ the ra- 
tio cto/cti increases monotonically with the filtering scale 
such that, even in the limit Rs — > 00 {y — \ 00) where 
the spatial bias is linear (b^(r) rs &„), the peak velocities 
remain fc-dependent. In general, the linear bias approx- 
imation 6n pk = b u 5s with unbiased velocities v pk = v s 
will provide a good description of the large-scale proper- 
ties of density peaks only when k <C min[y / 6^/5^, o^/do]. 
For density peaks of height v > 1, the square root ap- 
proximately is ((7i/(To)(^/v / 3)- The above condition thus 
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FIG. 1: Bias parameters <rob v , o"2&c an d ratio of spectral mo- 
ments <ro/o"i as a function of the filtering scale. Results are 
shown for density maxima of height u = S sc /ao at redshift 
z — 0.5. Dotted curves denote negative values. The linear 
spatial bias of peaks becomes scale independent in the limit 
v — > oo. However, the fc-dependence of the velocity bias, 
which is controlled by ao/ai, remains and even increases with 
the peak height. 



becomes k <C a\ja^. For the density maxima consid- 
ered here, this implies that the fc-independent linear bias 
approximation will be accurate for fc <C 0.1 hMpc^ 1 . 

The three-dimensional velocity dispersion of peaks is 
known to be smaller than that of the dark matter 



gradient of the density vanishes (by definition), and so 
v p k{xpk) — v(x p ) l ). The peak velocity dispersion is lower 
than that of the mass because large scale flows are more 
likely to be directed towards peaks than to be oriented 
randomly. This illustrates an important point: peaks are 
biased tracers which move with the dark matter flows - so 
although there is no physical bias in the velocities, there 
is a statistical bias which arises from the spatial bias. 
In the case of peaks, the spatial bias implies that 6 V ei 
is fc-dependent, and this introduces fc-dependence into a 
number of peak- velocity statistics (we provide an explicit 
calculation of this in equation I3"f51 below) . This is almost 
certainly true in general: fc-dependence of the spatial bias 
will lead to fc-dependence of 6 V ci even if the tracers flow 
with the dark matter. Note, however, that this is not 
a necessary condition since, for the highest peaks, the 
velocity bias remains scale dependent even though the 
spatial bias has no fc-dependence. 

The equality (6 V ci) = (fr V ci) does not uniquely constrain 
the velocity bias. For instance, the choice b vc \(k) = 
1 — (<J 2 _ 1 / 2 / cr o) k a l so nas (Kci) — (^vci)- However, if 
we think of the velocity bias as a real space operator 
b vc \(x) that maps a vector (velocity) field onto another 
vector field, then for homogeneous and isotropic random 
fields b VQ \(x) must transform as a scalar under rotations. 
Hence, it must be built from powers of the Laplacian V 2 , 
and this brings down a factor of k 2 upon a Fourier trans- 
formation. Therefore, we generically expect the lowest 
order fc-dependence to scale as 6 ve i(fc) = 1 — R^, c \k 2 (for 
some constant i? V oi), at least for tracers whose spatial 
bias relation can be expressed as a local mapping of the 
(smoothed) density and its derivatives. 

Before concluding, we emphasize that 6 p k(fc) and 
6 vc i(fc) are first order bias parameters. We expect con- 
tributions from higher order spatial and velocity bias pa- 
rameters to become more important as k increases, but 
calculating them is beyond the scope of this paper. 



vpk 



= ^ 1 (l-7o 2 ) 



(22) 



Notice that Eq. ([21]) for the peak velocity bias yields the 
same number, 



vpk 



1 

2^2 



dkP Ss (k)b 2 vcl (k), 



as it should, but that 



(23) 



(24) 



also! If we regard the integral over one power of & V ci, say 
(& vc i), as the peak-dark matter velocity variance at the 
same point (when smoothed on the scale of the peak), 
then the fact that (6 ve i) = (fr vc i) indicates that, at the 
position of the peak, the velocities of the peak and the 
mass are the same. This can also be seen in the average 
bias relation, Eq. (QJ: at the position of the peak the 



III. REDSHIFT SPACE CLUSTERING OF 
DENSITY MAXIMA 

We derive three estimates of the redshift space clus- 
tering of peaks. The first generalizes the formulation of 
based on linear theory of gravitational instability; it 
furnishes a simple estimate of the power spectrum. The 
second extends the probabilistic interpretation of it 
provides an expression for the correlation function. Ref- 
erence [45j has emphasized that, within the context of 
linear theory, this description of the correlation function 
is exact whereas that of [2J is only approximate. Ana- 
lytic approximations based on a probabilistic treatment 
lead to terms which, upon Fourier transforming to ob- 
tain the power spectrum, are lacking in the approach of 
0. This has recently been emphasized by [l(|. Finally, 
our third estimate shows that if one Fourier transforms 
at an earlier stage in the analysis, one obtains a slightly 
more intuitive expression for the redshift space power. 
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We examine in detail the implications of this approach 
for density peaks, despite the fact that much of this was 
already done by [2(| , for the reasons stated in the Intro- 
duction. 



A. Simple estimate of redshift space clustering 

The redshift space coordinate (also in /i _1 Mpc since 
velocities are in unit of length) is given by s = (sm,Sj_), 



s = x + f [v(x) ■ z] z . 



(25) 



where z is the unit vector along the line of sight. There- 
fore, at the lowest order, the redshift space density con- 
trast is related to that in real space by 



S s {k, t i) = S(k) + ffi 2 e(k) 



(26) 



where fi is the cosine of the angle with the line of sight 
@. For peaks, this becomes 



Sn s pk (k, n) = <5n pk (fc) + f/j? pk (fc) 

= b pk (k)S(k) + ffb vcl (k)S(k) 
b vcl (k) 



1 



6pk(fc) 



Sn pk (k) 



(27) 



upon insertion of the peak bias relation (j4]). Note that 
//i 2 is now multiplied by a fc-dependent factor. 

Using the former relation, the calculation of the red- 
shift space power spectra at leading order is straightfor- 
ward and yields 



(6pk(fc) + [M*0 + M*01 ff 

+ & ve l(fc)/V) Ps(k) 

(K + b ( k 2 + b vcl (k) JV) 2 P S (k) 



(28) 
(29) 



(the reason for introducing the superscript will become 
clear shortly). However, the corresponding expressions 
for the redshift space correlations are lengthy; we provide 
them later in this Section. 

It is conventional to write the redshift space power 
spectrum in terms of the real-space one, 



fovcl(fc) 

b P k(k) 



ff 



Ppk(fc). (30) 



The window functions cancel out, leaving us with 



frvcl(fc) 

bpk(k) 



1-k 2 



-by 



l 



(fc< 1) 



(TO 



1 / V 
1 — — rr 



(31) 



k 2 \a bv 



(fc> 1) 



Hence, unless care is taken, this will lead to constraints 
which depend on k even in the limit v 3> 1 where b v — > 
v/ cra- 



Ti. Probabilistic treatment 

In linear theory, the redshift space two-point corre- 
lation function £ s is related to that in real space by a 
convolution of the two-point correlation function in real 
space, £(r), with the probability distribution for veloci- 
ties along the line of sight 0, Q : 



i+r(sn,s ± ) = 

where 



dy K(y) 



2tt fa 12 {r) 



exp 



(s\\ -y) 2 



K(y) = l + £(r) + (I) 



0-12 (r) \fai 2 {r 



IS 



(32) 



»(r) 



'12 



(r) 



( s \\-y 
VM2W 



(33) 



Here, v\2(r) and 012(7") are the mean and dispersion of 
the pairwise velocity distribution of pairs separated by r 



in real-space (note that 



y 2 +s\ 



As emphasized 



by [45| , within the context of linear theory and the plane- 
parallel approximation, this expression is exact. This for- 
mulation is usually referred to as the "streaming" model. 
It should be noted that random pairs in real space are 
mapped to real space differently at different separation 
r because the pairwise velocity distribution depends on 
scale 0. 

Equation ([521 can be generalized to give 1 + £p k , the 
redshift space correlation function of peaks, simply by 
replacing V12 and CT12 with the expressions appropriate 
for peaks 20]. At first order, these are 



Parameter constraints are then derived from the angular 
dependence of Pp k , under the assumption of linear scale 
independent bias, for which b^ = and & V oi = 1, so the 
term which multiplies /i 2 is f /b u , and b v is assumed to 
be a constant. Our analysis shows that, for peaks, this 
prefactor is fc-dependent, and it depends on peak height. 



«i2(r,M) = [1 + f P k] 1 



2fc, ( ^[ 1/2) ~ t 1/2) 
a i 



+ 2b c ( /2) 



£1 



(1/2) 



Liifi) , (34) 



7 



2 / 
3^ 



2\ 2 



3^2 1^0 ^SQ 



j(-l) 



2C: 



(o) 



9 S2 



(35) 



where /i = r • z is the cosine of the angle between the 
line of separation and the line of sight, and the Li(/i) are 
Legendre Polynomials [4o| . 

Appendix El demonstrates that Eq. exactly repro- 
duces the result of a lengthy derivation based on the peak 
constraint. Note however, that it can be derived simply 
from setting 



_ ((1 + Sn pk:1 )(l + 5n p k,2)Opk,i - v p m) ' z ) 



((1 + 5n pk ,i)(l + <5w P k,2)) 



(36) 

where the subscripts 1 and 2 indicate positions separated 
by r, and the average is over all peak pairs with sepa- 
ration r. The correspondence with equation Q34[) can be 
seen by noting that in fc-space, the spatial bias from <5n p k 
is the sum of two terms, one of which is proportional to 
k 2 (equation [7]) and the velocity bias (equation [2Tj) in- 
troduces additional k 2 terms which come with factors of 
(cto/ci) 2 . Each additional factor of k 2 changes £^ to 

An+l) 

si 

The first term on the right-hand side of Eq. (l3"5j) is 
twice the (one-dimensional) velocity dispersion of peaks; 
recall that it is reduced by a factor of 1 — -f 2 relative to 
that of the dark matter (see Eq. 



1. Approximating the integral 



where all quantities in the right hand side are evaluated 
at s and primes denote derivatives with respect to sn (re- 



call that s = 



Eq. (|37j) describes the large-scale 



limit of the redshift space correlation function, in which 
derivatives of the real space correlation and pairwise mo- 
ments (i.e. the distortions) are small [l(|. When applied 
to dark matter rather than density peaks, the Fourier 
transform of the first three terms on the right-hand side 
yields Kaiser's formula 0, E^|. The fourth term arises 
because the pairwise velocity dispersion does not vanish 
even in the large scale limit [10(. We show below that 
Fourier transforming the analogous terms for peaks gives 
Eq. ([291). 



The derivatives of £ P k(s), vi 2 (s) and cr 2 2 {s) with re- 
spect to the line of sight distance sn can be evaluated us- 
ing ds/ds\\ = [i and d[ijds\\ = (1 — fj,)s , which follow 
from the fact that s 2 — s 2j r s\. The following relations 
are useful: 



dsn 



ds 2 

sf^iM 



= 4sf +1/2) ^) + k" +1/2) 



ds 2 



y/r 1] U»)-y/2 n+1) L2(,) 



+ 15 4o 



(n+l) 



(38) 



When a 12 < S||, then the Gaussian term in the ex- As a rulGj terms in An) appear always multiplied by 
pression above will be sharply peaked around y = s\\. the Legendre polynomial of order I. The lowest even 
Expanding f, v 12 and a 12 about their redshift space val- polynomials are L (aO = 1, L 2 (jjl) = (3^ 2 - l)/2 and 
ues yields L ^ = (35^ 4 -30m 2 + 3)/8. The calculation of the red- 

shift space correlation of peaks £ p k( s iM) is now straight- 
, (37) forward. Adding all terms together, we find 
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— rso 
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(39) 



As can be seen, there are harmonics up to ^''(•s) which arise from the second derivative !;"(s) in Eq. ([37]). These 
terms are significant only at distances less than a few smoothing radii and across the baryon acoustic feature where 
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the density correlation ^q -* changes rapidly [29j . Furthermore, terms linear in / arise only from the derivative of the 
pairwise velocity, — fv[ 2 (s). 

The redshift space power spectrum Pp k (fc,/z) in this approximation is obtained simply by Fourier transforming 
Eq. (|39|) . For the sake of completeness, 
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^B(k) + A -B\k)^ 2 -fk-{l- 1 l)a\ 



l + ^(fc) + ^ 2 (fc)-i/ 2 fc 2 (l- 7o 2 )^ 1 



Ppk(fe) 



1 k2 -ff 2 (l-i 2 )^ lPp]l (k), 



1 + B(k)fi 2 P pk (fc) 



(40) 



where the linear redshift distortion parameter 

Aei(fc) 



B(k) = f 



bpk(k) 



(41) 



is scale dependent. Recall that & p k(fc) and b ve \(k) were 
defined in equations and (|21[) . Thus, except for the 
second term in the last equality, the above result exactly 
matches our simple estimate, Eq. (|2^|) . 

Notice especially that, for linearly biased tracers, red- 
shift space distortions are used to estimate j3 = f/b. The 
analogous quantity for peaks, B(k), is fc-dependent. We 
will consider the implications of this in the next section. 



2. A different approximation 

A more intuitive approximation to the exact result that 
is reached upon performing the integral in Eq. (|32p can 
be obtained by Fourier transforming it in the first place. 
We write 

exp (— ifc • s) = exp (— ifcj^ • sj_) exp (— ik^y^ 

x exp[-ifc|| (s\\-y)] (42) 
= exp (— ik ■ r) exp [ — ifc|| (s|| — y)] , 

and then rearrange the order of the integrals so that the 
integration over sii — y is done first. Next, we use the fact 
that 



dte -t 2 /2 e -ifc« 



-k 2 /2 



(43) 



dtte-* /2 e" ifct =-ifce- fc ' 2 

Jdtt 2 e- t2 / 2 e- ikt = (l-fc 2 )e- fc2 / 2 , 

to express the result of the integral over s 
y as exp[— f 2 k 2 a 2 2 {r) /2] times other factors, 
nally, we recast this term as exp[— / 2 fc 2 cr 2 2 (oo)/2] 
times exp[— f 2 k 2 (a 2 2 {r) — <r 2 2 (oo))/2] which for small 
fc|| is approximately exp[— / 2 fc 2 cr 2 2 (oo)/2] x [1 — 
/ 2 fc 2 (er 2 2 (r) — er 2 2 (oo))/2]. Thus, we generically expect 



Fi- 



the redshift space power spectrum to take the form 
cxp[— / 2 fc 2 CT 2 2 (oo)/2] times other factors. A little algebra 
shows that, to leading order, these factors are precisely 
those given by equation ([29)) . giving 



p pk( k ^) = ex P 



-fk 2 a 2 vclf i 2 



(44) 



Here, P£® is given by Eq. (f30"T) . We have also used the 
fact that, except in pathological cases, the pairwise dis- 
persion at very large separation is simply twice the one- 
dimensional velocity dispersion of single particles, er ve i 
(= 0Vpk/3 f° r peaks), in units of aHf. Our notation is 
purposely kept general to emphasize that these results 
apply to any tracers of the linear density field. 

Our equation (|4"4")) corrects a number of important er- 
rors in previous analyses [2(1 S3] ■ In addition, expand- 
ing the Gaussian smoothing term shows the origin of the 
extra terms identified in the previous subsection (those 
highlighted by [l0(). Finally, the form of our expression 
reflects the fact that the associated correlation function 
can be written as a convolution of the original expres- 
sion (which is the Fourier transform of Eq. [29)) with 
a Gaussian in the line of sight direction: 



+ 00 



ds\\ G 



5||,CTl 2 (oo) 



£pk(S-L,S|| 



*\\) ■ 

(45) 

Now the meaning of our notation should be clear: the 
superscript refers to the limit in which the dispersion 
of the Gaussian smoothing term is vanishingly small. We 
note that this form for £ s was shown to be appropriate for 
the dark matter, without using any Fourier-space analy- 
sis [45[ ; our analysis demonstrates that it carries through 
for peaks as well. The interesting subtlety brought by 
density peaks is that the amplitude of the damping term 
CT vo i is related to the form of & V oi (Eq. |2U| . 

The functional form of our equation (|4"4l has been 
studied previously in the context of modelling nonlinear 
corrections to the redshift space power of linearly biased 
tracers Q, although there the assumption was that 6 p k 
is constant and & V oi is unity. We will discuss the effects 
of nonlinearities shortly. For completeness here, we sim- 
ply borrow all that previous analysis to show the effect 
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that the Gaussian smoothing has on the Fourier space 
multipoles. These can be written as 

j^L = A ( K ) + ~A x {K)B{k) + -A 2 {n)B 2 {k) , (46) 



Vj(k) = 5 
Ppk(fc) 2 

+ 

and 

Vl(k) _ 63 



[A^k) - A (k)} 



3A 2 (k) - -Ai(k) 



— A 3 (k) - -A 2 (k) 



B 2 (k) , 



B{k) 
(47) 



Ppk(fe) 



A 2 (k) - -jAi(k) + y a o(k) 



+ 



(48) 

B(k) 
B 2 (k) , 



where k = fko~ vc \ and the coefficients Ai(k) are recur- 
sively defined as 



A (k) 



Hi erf(/t) 
2 ~ 
(2£+l) 



(49) 



(2£+l) 



The final approximations assume k <C 1. Thus, the low- 
est order corrections to P^pk/Ppk and ^pk/^pk are pro- 
portional to — k 2 . 



C. Nonlinear evolution 

There are four reasons why nonlinear evolution will 
act to change the expressions above [l(| HI] : one is re- 
lated to the change in the bias parameters, and the three 
others have to do with the effect of peculiar velocities. 
Gravitational motions are expected to relate a scale in- 
dependent, deterministic linear bias parameter in the ini- 
tial (Lagrangian) field to the evolved (Eulerian) bias ac- 
cording to b^ ul = 1 + b v [12] ■ Ignoring the fact that 
might also evolve, we follow common practice and as- 
sume = 1 + b v + b^k 2 , even though we suspect that 
bp^ 1 = b vc \ + b v + b^k 2 (This issue will be thoroughly ex- 
plored in a forthcoming paper) . Note that we have omit- 
ted the smoothing window for brevity, but it effectively 
makes little difference at scales A: -1 ^> Rs- Regarding 
the peak motions, we first assume that virial velocities 
within peaks or halos will increase 1712(00); these are re- 
sponsible for the fingers-of-god [54j, l55| seen in galaxy 
surveys. Secondly, halo/peak motions may not closely 
follow linear theory, but this effect is expected to be less 
dramatic (56| . Finally, the real space power spectrum will 
also be modified as a result of the linear theory motions 

MM- 



For reasons we describe below, nonlinear effects may 
be approximated by setting 



p p s k (fc, m) = P$(k, m) v^k, V 2 ) v vil (k, m 2 ) 



(50) 



where P*" is given by Eq. (|30[) with 6 p k replaced by its 
nonlinear version. The filters V q \ and V v i r are supposed 
to reflect the quasi-linear and virial corrections to the 
non-damped linear theory expression, respectively. The 
exact functional form of V q \(k, /i 2 ) depends upon the dis- 
tribution of pairwise velocities. However, motivated by 
results from perturbation theory [57l l6ll |62| , we set 



V q i(k,Li 2 ) = cxp 



-k 



.2/1 2\ 7 2 2 



(1 + /)V 



(51) 

This takes into account both the smearing of linear power 
caused by linear theory displacements and the damping 
due to the linear pairwise velocity dispersion. In princi- 
ple, the reduction in linear power should be somewhat 
mitigated by the addition of nonlinear mode-coupling 
terms of the sort discussed by 60]. However, we will 
ignore these terms in what follows. The last multiplica- 
tive factor V v i T accounts for the damping of redshift space 
power due to nonlinear virial motions within halos (as- 
sumed uncorrelated with the large-scale flows). If the 
mass range is small (i.e. if the peaks cover a small range 
in v), then V v - a = exp(— k 2 ^i 2 a 2 ir ), where <7 v i r depends 
on the halo or peak mass, should be a good approxima- 
tion. If the mass range is broad, then an exponential 
distribution may be more appropriate (53|, leading to 
V n \ = [1 + k 2 a 2 r{ At 2 l~ 2 . Removing fingers-of-god from 
a survey [e.g. |j5l| is equivalent to setting cr V i r — > or 
Kir -t 1. 

Here and henceforth, we will assume that V v i r is a 
Gaussian smoothing kernel. This implies that the Fourier 
space multipoles Vf{k) are given by Eqs (|4"6"|) . (|4"7| and 
with 



(52) 



n^k^a 2 vch f(2 + f) + l 
upon making the replacement 

V!(k) i pf(fc) 



P P k(fc) P pk (fc)e- 



fc 2 cr 2 



(53) 



on the left-hand side. We will now illustrate the effect 
of the biasing relation Eq. (j4|) on the 2-point correlation 
through a comparison between density peaks and linearly 
biased tracers. 



D. Comparison between density peaks and linearly 
biased tracers 

For linearly biased tracers, b^ = 0, 70 = and all 
the terms involving ctq / ci vanish. The pairwise statistics 
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FIG. 2: A comparison between the redshift space multipoles of the correlation function of density maxima and linearly biased 
tracers. (Dotted lines denote negative values.) The peaks were identified in the density field when smoothed with a Gaussian 
filter of characteristic scale Rs = 2.5 (left panel) and 4 /i _1 Mpc (right panel). This corresponds to a mass scale Ms = 1.9 x 10 13 
and 7.8 x 10 13 Mq/Zj, respectively. The associated peak height and bias parameters quoted in each panel assume a redshift 
z = 0.5. The linear biased tracers are required to have the same value of 6^ ul . The peak biasing relation enhances the monopole 
and the quadrupole around the BAO scale relative to that of linearly biased tracers, and induces significant scale dependence 
in the hexadecapole at s < 100 Mpc. 



simplify to 

-2b^[- 1/2 \r) 
vi2(r,H) = — - ■ - , , — Ltfjl) , 



0"12 (»*>M) = 



So 



(-i) 



Mm) 



(54) 
(55) 



Setting j3 = f/b„, we recover the linear theory prediction 
of [H plus a contribution from the large-scale limit of of 2 
(which underestimates the true effect since we neglect 
nonlinear corrections to the velocity dispersion), 



(56) 



2 AO) 



L2M 



1 



f2^2 .2 C (1) 



Note that Eq. ([56)1 implicitly assumes that the peculiar 
velocities of linear tracers match locally that of the mat- 
ter. To account for the nonlinear evolution, we will also 
adopt the prescription b v — > b*j ul = l + b v . 

The explicit Legendre decomposition £ s (s, /i) = 
5j£!(s)X^Gu) °f the redshift space correlation function 
can be read off from Eqs. (|39l) and (I5"6"l) . For illustra- 
tion, the multipoles £f(s, v) are plotted in the left and 



right panel of Figf2] for density maxima identified at the 
smoothing scale Rs = 2.5 and 4 /i _1 Mpc, respectively. 
These functions are compared to those of linearly biased 
tracers with same value of &^ ul , namely b^ ul — 2.0 and 
3.8 respectively. It is important to note that, in the lat- 
ter case, the density field is not smoothed (in practice 
we use Rs =0.1 /iMpc -1 ). Furthermore, we have also 
neglected the contribution <j V i r from virialized motions 
to the velocity dispersion and assumed cr^ ol = cr^ pk /3 for 



.11 /i^Mpc) 2 , for 



the peaks and a 2 L 1 /3, with a 2 L 1 - 
the linearly biased tracers. 

As recognized in [29j, the nonlinear local biasing re- 
lation Eq. (j4} amplifies the contrast of the (real space) 
baryon acoustic signature of density maxima relative to 
that of linearly biased tracers. A similar enhancement 
is also observed in the baryonic acoustic signature of 
dark matter halos in very large cosmological simulations 
[HI, H|[ . As can be seen in FigfSl this amplification is also 
present in redshift space. In this case however, both the 
monopole and the quadrupole of ^ k are affected by the 
nonlinear peak biasing across the baryon acoustic oscil- 
lation (BAO). At distances s - 100 - 110 /i _1 Mpc, the 
quadrupole of the redshift space peak correlation £* k is 
indeed more negative, damping thereby the correlation 
in the radial direction (/1 w 1) and increasing it in the 
perpendicular direction (/i w 0). 

This is more clearly seen in Fig. [31 which compares 
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s / fr'Mpc s / fr'Mpc 

FIG. 3: Angular dependence of the redshift space correlation for the density maxima (solid curves) and linearly biased tracers 
(dashed curves) considered in Fig. [2] Results are shown for a separation vector oriented in the direction parallel (fi = 1) and 
transverse (fi — 0) to the line of sight. Perpendicular to the line of sight, the contrast of the acoustic peak is more pronounced 
in the correlation of density peak whereas, in the radial direction, it is comparable to that of linearly biased tracers. 



the redshift space correlation of density peaks and lin- 
ear tracers in the direction parallel and transverse to the 
line of sight axis. Relative to the baryon acoustic peak of 
linearly biased tracers, the BAO of density maxima is en- 
hanced in the direction perpendicular to the line of sight 
while somewhat distorted in the radial direction. The 
physical origin of this effect presumably is peak-peak ex- 
clusion. Namely, while as discussed in |29| the spatial 
bias of peaks enhances the contrast of the BAO in the 
real space correlation, peak-peak exclusion suppresses the 
infall of peak pairs onto the (slightly overdense) BAO 
shell at radius s ps 105 /i _1 Mpc. In redshift space, 
this amounts to a reduction of the BAO contrast along 
the line of sight. The dispersion term a^ipo)^" further 
smoothes the BAO and shifts the position of the local 
maximum in that direction, but leaves the baryon wiggle 
unchanged in the transverse direction. The amount of 
smoothing depends on the exact value of cr 2 el . Another 
striking feature of Fig. [3] is the strong suppression of the 
redshift space correlation and the sharpening the acous- 
tic peak along the line of sight due to linear coherent 
infall 0,1. 

On scales less than the BAO ring, the contribution of 
the pairwise velocity dispersion increases with decreasing 
separation until it reverses the sign of the quadrupole 
at separation ~ 10 — 20 fe~ 1 Mpc and stretches struc- 
tures along the line of sight [l(|. Although the veloc- 
ity dispersion of the density peaks is smaller than that 
of the linear biased tracers, peak-peak exclusion makes 
the effect stronger. On those scales, the contribution 



of the term 2b^ ul b^Q + b^ffl becomes comparable to 

(6^ ul ) 2 £o°^ and steepens the profile of the angle-averaged 
correlation As a result, the monopole for the density 
peaks can be larger by a few tens of per cent at separa- 
tion s < 10 /i~ 1 Mpc relative to that of linearly biased 
tracers. Note that small-scale halo exclusion is not prop- 
erly accounted for in our treatment since we consider £ p k 
at first order only. Nevertheless, we expect that, while in 
real space peak-peak exclusion leads to a deficit of pairs 
at distance s < Rs, in redshift space the suppression 
may be weaker because peaks tend to move toward each 
other. 

Fig. @] displays the Fourier space multipoles Vi(k) in 
unit of -Ppk(fc) for the peaks and linear tracers considered 
above. To emphasize the importance of the exponen- 
tial damping, results are shown with and without the 
smearing caused by quasi-linear and virialized motions. 
While for the linearly biased tracers the distortion pa- 
rameter B(k) = f/b^ ul is a constant, for peaks B(k) is 
/c-dependent and, therefore, induces a scale dependence 
in the multipoles even when the pairwise velocity disper- 
sion is negligible. In this limit (a ve \ = 0), for the lin- 
early biased tracers the ratios Vt{k)/ P p k(k) are constant 
(as in the original Kaiser formula), whereas for peaks 
they decay rapidly to reach 1 + 2£?(oo)/3 + £> 2 (oo)/5, 
4£(oo)/3 + 4£ 2 (oo)/7 and 82? 2 (oo)/35 when I = 0, 2 and 
4, respectively. Here, B(oo) = — {v/gq — b® )^ is the 
value of B(k) in the limit k — > oo (see Eq. l3"Tj) . The 
difference between peaks and linear tracers is largest in 
the hexadecapole and increases with mass scale. For in- 
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k / hMpc" 1 k / hMpc" 1 

FIG. 4: Fourier space multipoles "Pf in unit of P pk (k) as a function of wavenumber for the density peaks and for the linearly 
biased tracers with same value of b^ ul . The dotted-dashed (peaks) and dotted curves (linear tracers) are the results without 
including the velocity damping kernel (assumed to be a Gaussian, see Sec. IIIIC)l . while the solid (peaks) and dashed (linear 
tracers) curves represent the multipoles when the Gaussian dispersion is included. For density peaks, V£(k) exhibit a strong 
fc-dependence even upon removal of the damping term. 



stance, the fractional deviation is 5 per cent at wavenum- 
ber k ~ 0.037 /iMpc -1 and « 0.027 /iMpc" 1 for the peaks 
identified at filtering scale R$ — 2.5 and 4 /j _1 Mpc, re- 
spectively. 

When the Gaussian damping term, Eq. (j5"Tj) , is in- 
cluded, the behaviour of the Fourier space multipoles of 
density peaks (solid curves) and linear tracers (dashed 
curves) becomes similar at small-scale: they damp to 
zero like the coefficients A^k) defined in Eq. f4"9"]) . Still, 
significant deviations persist on scale k > 0.01 /iMpc -1 
due to the fc-dependence of B(k) and unequal velocity 
dispersions. 



IV. COSMOLOGICAL IMPLICATIONS 

A. Estimating the growth rate / 

Following @, E3, HE] i the redshift space power spectrum 
of density peaks can also be written as 



P" k (fc) = P pk (fc) + 2/, 2 P pMpk (fc) + //P 0pk (AO F{k, m 2 ) 



dependent of the spatial bias and directly measures the 
matter velocity power spectrum provided there is no ve- 
locity bias. Owing to the angular dependence, a mea- 
surement of the velocity power spectrum furnishes an es- 
timate of the linear growth rate fag oc dD/dXna that is 
not affected by the spatial bias. 

Although the hexadecapole does not depend upon the 
spatial bias, it may be noisier than the monopole and 
dipole, so this has motivated the search for other com- 
binations of Prj and Pi which may be more robust [8|. 
Reference [3l[ showed that Pq and P| can be used to de- 
rive an estimate of the velocity power spectrum Pe k (fe) 
when the density and velocity fields are perfectly corre- 
lated, namely, when the cross-correlation coefficient 



p2 
pk,i 



Ppk(fc)PvW 



(58) 



is unity. For example, when smoothing is ignored, then 



P : 



245 
l8~ J 




(59) 



(57) [31| is proportional to f Pg(k) when velocities are un- 



where P p k,e pk (fc) and Pe pk (fc) are the peak-velocity biased. Here, P| = P 2 S /P S - For peaks rg(fc) = 1 in- 
and velocity- velocity power spectra. Here, P p k, 
Ppk,6» pk and Pe pk are linear spectra and P(fc,/x 2 ) = 
V q i(fc,/z 2 )Kir(fc,^ 2 ) describes both the quasi-linear 
damping and the smearing from the small-scale velocity 
dispersion. As noted in [3l|, [H, [6(1 [67[, Pe pk {k) is in- 



deed holds at the lowest order, even though the linear 
spatial and velocity bias b p ^(k) and b ve i(k) are scale de- 
pendent. However, the velocity power spectrum now 
is Pe pk {k) = f 2 bl cl (k)Ps(k), so there is an extra k- 
dependence associated with the estimator P. Since this 
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could be interpreted erroneously as a signature of modi- 
fied dark energy or gravity, any scale dependent velocity 
bias (a scale independent bias may also be present if the 
tracers do not move with the matter) will limit the in- 
formation that can be recovered about the growth factor 
[3ll [52l . [66|. For peaks, the velocity bias is b ve i(k) < 1, 
and it converges towards unity (i.e. unbiased velocities) 
in the limit k — > 0. At the first order, the deviation from 
unity is controlled by ao/ci (Eq. [21]) so that, at fixed 
wavenumber, b ve i(k) decreases with increasing mass scale 
(see Fig. □ For M s = 1.9 and 7.8 x 10 13 M & /h consid- 
ered here, Pg pk (k) is suppressed by « 5 and 9 per cent 
at wavenumber k ~ 0.05 /iMpc -1 , respectively. The pre- 
dicted fc-dependence is smaller than current constraints 
on the growth rate [52, 68]. Furthermore, numerical sim- 
ulations to date show that the power spectrum of dark 
matter halo velocities is consistent with f 2 Ps(k) within 
10 per cent at wavenumber k < 0.1 /iMpc -1 [3l|. Never- 
theless, since forthcoming large-scale galaxy surveys will 
dramatically improve constraints on the growth factor 
(down to the percent level), it is interesting to assess the 
extent to which a fc-dependent bias would degrade the 
constraint on the growth rate. 



B. Error forecast with a fc-dependent velocity bias 

To this purpose, we use the Fisher based formalism 
developed in [66|]. For Gaussian random fields, the Fisher 
matrix for a set of parameters {pi} is [TOL l71j 



d 3 k f d\nP \ ( d\nP \ 
(2^ \®w) {d^p~J M ) 



(60) 



where P is the power spectrum and individual wavemode 
contributions are weighted by the effective volume [69] 



V eS (k) = V 



nP 



1 + nP 



(61) 



which depends upon the surveyed volume V and the 
number density n of the tracers (assumed homogeneously 
distributed). To illustrate, we assume the linear, plane- 
parallel approximation and consider the model 



P s (k,fi)= b pk (k) + fb vcl (k)fi 2 P s (k), 



(62) 



where b pk (k) 



b^k 2 (we drop the superscript Eul 



for brevity), and b ve i(k) = 1 — R 2 c \k 2 (for some i? vo i) 
are motivated by the functional form of the spatial and 
velocity bias of density peaks (c.f., Section lll Cl andl ll Fl) . 

In what follows, we fix the shape and amplitude of 
the matter power spectrum (i.e. the fractional error 
on fag is equal to that on /) and consider the four- 
parameter set {fry, i? vc i, /}. Our fiducial model has 
(b vi , i?voh /) = (1,16,3,0.46). The values of b^ and 
i? V ci closely correspond to those of density peaks identi- 
fied at the mass scale 1.9 x 10 13 M Q //i. Derivatives of 



the logarithm of the power with respect to the parame- 
ters are computed easily: 



<91nP 



<91nP 



2k 2 



db v (b pk + fb vc ifi 2 )' db c (6 pk + fb vcl n 2 ) ' 
<91nP -4/i? vclM 2 fc 2 



<9i?vel (6pk + fKclH 2 ) 

d\nP _ 2b ve m 2 
df {bpk + /fovoiM 2 ) ' 



(63) 



We integrate over wavenumbers from fc m ; n ~ tt/V 1 ^, 
where V is the volume of the survey, up to a maximum 
wavenumber k = 0.1 /iMpc , above which nonlinear 
effects are expected to become important (66|. 

In order to illustrate the effect of including a fc- 
dependent velocity bias into the analysis, we initially 
set 6 V ei = 1 (i.e. ignore i? V ci) and compute the Fisher 
matrix for b„, b^ and / solely. For a survey of volume 
V = 10 /i~ 3 Gpc 3 at redshift z = 0, we find a fractional 
marginalized error of Sf/f =1.6% in the limit nP 1 
of negligible shot noise (In practice, a suitable weighting 
of galaxies may help approaching this limit [73[ ). For a 
number density n = 5x 10 -4 and 10 -4 h 3 Mpc~ 3 , the con- 
straint weakens to 1.9% and 2.9%, respectively. (These 
values are consistent with those of [13].) Unsurprisingly, 
6„ and b^ are strongly anti-correlated (the correlation 
coefficient is r rs —0.8) because an increase in b v can be 
mostly compensated by a decrease in b^. However, while 
the correlation between 6„ and / is moderate (r < —0.5), 
b^ and / are weakly degenerate (r < —0.05). In other 
words, including a fc-dependent bias component b^k 2 has 
little effect on the uncertainty on /. Extending fc max be- 
yond 0.1 /iMpc -1 (where the shot noise becomes again 
important) can reduce the uncertainty on / (because the 

fractional error scales as fcmax 2 ), but this is at the price of 
having to model the smearing due to quasi-linear motions 
and small-scale velocities. 

Introducing the parameter i? vc i substantially increases 
the uncertainty on /. For the volume V and the aver- 
age number densities n considered above, the fractional 
marginalized uncertainty on the growth rate becomes 
Sf/f=A%, 4.4% and 6%, respectively. This can be traced 
to the strong correlation (r sw 0.9) between i? ve i and /. 
The error degradation reflects the fact that we are adding 
more freedom to the model. It does not depend upon the 
exact value of i? V ei- 

Are the constraints on f obtained using the multi- 
tracer method proposed in [751 ] affected in a similar way ? 
Reference [76j pointed out that several populations of 
differently biased tracers can achieve a much better de- 
termination of the growth rate than a single sample of 
objects. When power spectra are measured, calculating 
the Fisher matrix for multiple tracers requires summing 
over the distinct components of the inverse covariance 
matrix CJ S , where A, B label a different pair of tracer 
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populations, 



F; 



V 



A.B 



d 3 k (8P A \ ! 



AB 



3P B 

9 Pj 



(64) 



The calculation of the covariance matrix is straightfor- 
ward if one assumes that the noise term can be treated 
as an uncorrelated normal variate (66l . . For complete- 
ness, the diagonal and off-diagonal components of the 
covariance matrix are 16 61 



(C a 



and 



= 2P 2 a A 2 , (C abab ) = Pl h + P aa P bb N a N b (65) 



{C ab cd) — 2PabPcdj (C a abc) — ^PabPac 



(Cabac) 



PnhPn 



PanPhrN n 



(66) 



(C a 



PabPaa-N a , 



(C a a 



bb 



ZPZb 



where Pjj is an auto- or cross-power spectrum and N a 
I + l/(nP aa )]. We will consider the simplest case of two 
types of tracers, since the gains saturate rapidly as the 
number of samples increase 66] . In this case there are 3 
distinct measured power spectra 



P 



where a, & = 1,2 and the bias factors are = b„ + 

b^k 2 , b^J = 1 — (R^J) 2 ^. The power spectra are thus 
described by 7 parameters, and their derivatives are 
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Here, 5^ is the Kronecker delta. 

Fig. [5] shows the fractional marginalized error on / 
obtained by combining two different biased sample of 
the same survey volume. The constraints are shown as 
a function of the abundance of the second tracers with 
and without including a fc-dependent velocity bias (solid 
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FIG. 5: The fractional marginalized error 5f/f for a sur- 
vey volume V = 10 /i _3 Gpc 3 at z = obtained with two 
tracer populations : a high density, unbiased sample with 
ni = 10~ 2 /i 3 Mpc~ 3 and fc^ 1 ' = 1 and a second popu- 

(2) 

lation with varying number density ni and bias b v . We 
choose &1 2 ' = 1.4, 2 and 4 (curves from top to bottom). 
The constraints are shown as a function of ni assuming 
6^1 = b^j = 1 (dotted curves) and a fc-dependent velocity 
bias with = R^J = 3 /i -1 Mpc (solid curves). 



(2) 

and dotted curves, respectively). We set b v = 1.4, 2 
and 4 to facilitate the comparison with Fig. 3 of [6/ 
Although we have assumed the fiducial values R[ 

3 ft Mpc for simplicity, we may expect from 



(i) 



R 



(2) 
Vel 



the analysis done in the previous Section that i? V ei has 
some mass or bias dependence. E.g., R^J > R^J when 
b^ 3> bv . The marginalized error on / is, however, 
weakly dependent on the fiducial value of Note 
the considerable improvement in the constraint on / [in 
agreement with the findings of [66|, [76| . The smallest error 
is achieved with a large number density n-i and large rela- 
tive bias b^ /b^}\ (We have used values of b^ to simplify 
comparison with |66|.) However, including a fc-dependent 
velocity bias degrades the uncertainty on the growth rate 
roughly by a factor of two when n-z > 10~ 2 h 3 Mpc~ 3 , like 
in the single tracer case. The error degradation becomes 
increasingly severe as one goes to lower number densities. 

Although these constraints are only indicative (We 
have ignored the influence of cosmological parameters 
on the constraint [z3), our analysis demonstrates that 
allowing for a fc-dependent velocity bias (with the spe- 
cific functional form predicted by the peak model) has a 
large impact on the determination of the growth factor 
and, therefore, may possibly hamper our ability to distin- 
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guish between different dark energy or gravity scenarios 
[771 l78j . Therefore, despite the lack of current evidence 
for a fc-dependent velocity bias [3ljj , it seems prudent to 
study this possibility further with large cosmological sim- 
ulations. We hope the peak model can serve as a useful 
baseline with which to compare the simulations. 



V. STOCHASTICITY 
A. Cross-correlation coefficient 

The biasing eq. (j4j derived from the large-scale proper- 
ties of peak correlation functions is a mean bias relation 
that does not contain any information about stochastic- 
ity. Therefore, it is unsurprisingly deterministic like the 
local bias model considered by [39(, the main difference 
residing in the fact that peak biasing involves derivatives 
of the density field. Still, because of the discrete nature of 
density peaks, one can expect that the peak overdensity 
<5rt p k at location x generally be a random function of the 
underlying matter density (and its derivatives) in some 
neighbourhood of that point. We note that stochastic 
models of the form (5n p k(a:) = X[5s(x)] have been stud- 
ied in [79l - [8l| for instance. 

Computing the probability of X given 5 etc. is be- 
yond the scope of this paper (because it requires the full 
hierarchy of correlation functions). Still, it is instruc- 
tive to compute the cross-correlation coefficient to gain 
further understanding of the peak biasing model. The 
cross-correlation coefficient is defined as 



r 2 M 



Ppk(k)P s (k)' 



■(r) 



pk,S 



(r) 



£pk(r)fr(r) 



(72) 



in Fourier and configuration space, respectively. Ignor- 
ing the damping term, we find r c (k) = 1 for peaks even 
though the ratio Pp^^/Ps depends on k. Thus, a k- 
dependent bias at the linear order does not yield stochas- 
ticity in Fourier space. On the other hand, 



?(r) 



h 



(0) 



bi + 2b v b c & ) /$ ) + b 2 c &>/$ 



(2) /AO) 



(73) 



Therefore, although the bias is deterministic in Fourier 
space, it is generally stochastic and scale dependent in 
configuration space. However, when v 1 then — » 0, 
so r£ — > 1. Namely, in the high peak limit, the bias be- 
comes linear and deterministic in both Fourier and con- 
figuration space. 

The real space cross-correlation coefficient is shown 
in Fig. H as a function of comoving separation for the 
peaks identified at the smoothing radius Rs — 2.5 and 
4 /i _1 Mpc. is very close to unity at all separa- 
tions larger than a few smoothing radii, except around 
the baryonic bump and the zero crossing of the corre- 
lation function where it can be noticeably larger than 
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FIG. 6: The cross-correlation coefficient r%(r) for density 
peaks identified at the smoothing scale Rs = 2.5 and 
4 /i _1 Mpc (upper and lower panel, respectively). There is 
significant stochasticity only at the zero-crossings of the auto 
and cross-correlation functions and across the BAO, where 
the derivatives of the density correlation ^g 1 ' and £q are not 
negligible [IJ. 



unity. These findings seem to run contrary to the com- 
mon knowledge that \r^\ < 1. However, at separation 

r > 120 /i~ 1 Mpc, the fact that the zero-crossings of ^q™" 1 
do not generally coincide unavoidably implies \r^\ > 1, at 
least over some range of scales. Furthermore, at distance 
r ~ 90 /i _1 Mpc, the large values of are most plausi- 
bly traced to the baryon acoustic feature, which induces 

large oscillations in ^q 1 -* and f^ 2 -* across the BAO scale 
w 105 /i -1 Mpc (see Fig.l of [29|). At the level of a bias 
relation <5n p k = X[5s, V 2 5s,- ■■], this suggests that the 
scatter is strongly sensitive to V 2 ^. 

Figure [7] explores the behaviour of the cross-correlation 
coefficient when the underlying power spectrum is a fea- 
tureless power law spectrum, Ps{k) oc k ns . Results 
are presented as a function of the spectral index n s 
for a single value of the separation, r = 100 /i~ 1 Mpc. 
At fixed value of n s , the stochasticity is unsurprisingly 
larger for the relatively sparser peaks identified at scale 
Rs = 4 h~ 1 Mpc. Most importantly, the amount of 
stochasticity depends sensitively upon the shape of the 
matter power spectrum. Overall, decreases with larger 
values of the powerlaw exponent n s because the stochas- 
ticity rises as the relative amount of small-scale power 
increases. As can also be seen, is slightly larger than 
unity in the range —3 < n s < —2 and at the points of 
discontinuity n s = 0, 2 (which are marked as empty sym- 
bols). Although the effect is admittedly small and local- 
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ized in n s , this demonstrates that the cross-correlation 
coefficient can exceed unity also when the power spec- 
trum is scale-free. In Appendix |EJ we investigate the dis- 
continuities in more detail and provide quantitative esti- 
mates of the large-scale behaviour of the cross-correlation 
coefficient for a few values of n s . 



B. Evolution of stochastic bias 

As discussed in Sec lIII Cl gravitational evolution maps 
a scale independent, deterministic linear bias factor in the 
initial conditions onto a similar quantity in the evolved 
distribution [42[ . The scaling £>^ ul = l + b v also works for 
a fc-dependent deterministic bias. More precisely, 



bX(k,z)-l = b pk (k,z) = 



D(z)/D(z ) 



frpk(fc, Z ) 



D(z)/D(z ) 
(74) 

where D(z) is the linear theory growth factor [82|, |83| . 
This is easily understood if one recognizes that a peak 
of height b z S z , where S z is the linearly evolved field at z, 
could also have been written as having height bgdo where 
do is the field evolved to zq. The relation b z S z = boSo 
implies b = b z (S z /Sq) = b z D(z)/D(z ), from which the 
above expression is derived. Alternatively, notice that 
b oc l/a (z ) oc D(z)/D(z )/a (z) oc b z D(z)/ D(z ), 
so the factor D(z)/ D(z ) is simply converting from one 
choice of fiducial time to another. 

In configuration space, it has been argued that for lin- 
ear stochastic bias, 



bf u \z)r((z) 



bf»\z ) n (z )-l 



(75) 



D(z)/D(z ) 

[e.g. Hi)]. The corresponding expression for the evolution 



of bf u \r, z) itself is 



.Eul 



4 (z) 1 D 2 (z) = b^(z ) 2 D 2 (z ) + [D(z ) 
-2[D(z Q )-D(z)]D(z ) 



D{zf 



x bf ul < 



(76) 



This is a good model of bf for density peaks, pro- 
vided we interpret the denominator of Eq. ([73|) as 6|. 
Moreover, the numerator of this equation is similar 
to (the square of) an Eulerian bias factor minus one: 
([1 + b v + 6 C ^ 1) /^ 0) ] - l) 2 . This quantity clearly scales 
with the growth factor like its Fourier space analog. 
Hence, the real space evolution of the stochastic bias of 
density peaks is simple in spite of the additional scale 
dependence. Notice that both (k) and b® tend to 
unity at late times (even though they might effectively 
not reach this limit because the growth factors freeze out 
in ACDM-like models). In fact 6j does as well, but this 
is not as easy to see from our expressions. 
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FIG. 7: Cross-correlation coefficient r^(r) for powerlaw power 
spectra Ps(k) oc k" 3 as a function of the spectral index n s . 
Results are shown at a single separation r — 100 /i -1 Mpc, 
for density peaks identified at the smoothing scale Rs = 2.5 
(top) and 4 /i~ 1 Mpc (bottom). The insert is an enlarged view 
of in the range —3 < n s < — 1. Notice the discontinuities 
at n 3 — and 2, at which the cross-correlation coefficient is 
slightly larger than unity. 



C. Connection to previous work 

We mentioned earlier that peak bias and its evolution 
have been studied in simulations by [26|, [ijj . These au- 
thors found that the peak background split argument 
(Eq. [T5)) provides a good description of the large scale 
bias of the peaks extracted from their simulations. They 
also found that Eq. (|T4[) is in reasonable agreement with 
the evolution of this large scale bias. However, on smaller 
scales, the real space bias was found to be scale de- 
pendent and stochastic [4l[, two features which a peak- 
background split based analysis docs not model. Never- 
theless, Eqs f74"P and (|75|) were found to provide a good 
description of the evolution. The above analysis shows 
why. It would be interesting to see if our approach cor- 
rectly predicts the scale dependence of the bias. We defer 
this issue to a future work. 



VI. DISCUSSION AND CONCLUSIONS 

We have presented an extensive analysis of the Gaus- 
sian peak model. Density peaks are biased tracers of the 
underlying matter density field - on large scales this bias 
is scale independent - and we studied the limit in which 
this bias just starts to exhibit a scale dependence, both 
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in the spatial and velocity fields (Eq. [5]) . In almost all 
cases we presented a relatively straightforward analysis 
in the main text, which was sometimes backed up with 
detailed calculations in the Appendix. 

We showed that, in the large scale, scale independent 
limit, our expressions reduce to those of the peak back- 
ground split (Section III Dp . but in general the scale de- 
pendence in our model implies a much richer structure. 
For example, even though the peaks flow with the un- 
derlying field, their velocities appear to be biased (Sec- 
tion III F[) . In addition, we showed how this fc-dependent 
bias propagates into the analysis of redshift space distor- 
tions (Section IIIip . We derived an exact formula for the 
linear theory redshift space correlation function (Eq. 152]) , 
and then argued that it should be well- approximated by 
a simpler expression which has considerable intuitive ap- 
peal (Eq.HU). 

Our formula shows that redshift space distortions of 
peaks can be modelled i) using the same formula (and 
physics) as in reference 0], except that various terms 
now become fc-dependent, and ii) there is in addition 
a Gaussian smoothing term, which reflects the disper- 
sion of particle velocities in linear theory. Thus, our for- 
mula has the same form as the phenomenological rela- 
tion that is commonly used to model nonlinear effects, 
and which has been shown to provide increased accuracy 
when comparing theory with simulations. However, here, 
we demonstrated explicitly that this functional form is 
also part and parcel of linear theory. Kaiser's relation 
assumes that the smoothing term is unity (the fc <C 1 
limit) whereas Scoccimarro's formula [13], which is de- 
rived from the full Gaussian random field expression, is 
equivalent to expanding the Gaussian smoothing term 
and retaining only the monopole and quadrupole. Our 
result implies that linear theory can account for some 
of the effects that such a phenomenological model would 
otherwise ascribe to nonlinear evolution. 

We provided a crude treatment of nonlinear effects 
(Section IlllCj) . which, though not properly accounting 
for the nonlinear evolution of the matter density and 
velocity fields [S6j] nor for the mode-coupling contribu- 
tion induced by nonlinear gravitational clustering [13], 
illustrates how the new smoothing term, and the k- 
dependence of the (spatial and velocity) bias factors, im- 
pacts cosmological constraints from galaxy redshift sur- 
veys (Section livp. Our analysis showed that allowing for 
a fc-dependent velocity bias degrades constraints on the 
growth rate / by at least a factor of two. Large cosmolog- 
ical simulations will be needed to ascertain whether dark 
matter halos hosting the surveyed galaxies also exhibit 
a fc-dependent velocity bias. If they do, then improving 
the determination of / will lie in our ability to model this 
bias. 

We also used the peaks bias model to investigate the 
stochasticity of the bias and its evolution (Section [Vj) . 
We provided explicit expressions for the evolution of the 
scale dependent peaks bias and stochasticity, and argued 
that they helped to understand recent measurements of 



these quantities in numerical simulations. 

As regards the evolution of peak bias, it is interesting 
to consider the peak model in light of recent work on 
possible modification of gravity. In standard gravity, the 
linear theory growth factor is scale independent. There- 
fore, a peak retains its height when the initial density 
field is linearly evolved. However, in modified gravity 
models, the linear growth factor is fc-dependent. As a 
result, peaks in the initial field may not correspond to 
peaks in the linearly evolved field because the shape of 
the power spectra for the two Gaussian fields is different. 
This can also be seen directly by studying the (linear the- 
ory) motions of peaks. In such models, the bias of objects 
which coherently flow with the matter evolves just as it 
does in standard gravity 83]. Thus, whereas objects ini- 
tially placed at maxima of the density field will still move 
in accordance with the (modified) matter flows, gravita- 
tional motions will bring them to positions which are no 
longer local maxima. 

The question then arises as to whether it is the initial 
peaks, or those in the evolved field, which bear a closer 
resemblance to the galaxies and clusters we see today. 
Presumably it is the peaks which have managed to sur- 
vive from the initial time to the present which are the 
ones of most interest - the ones which are transients are 
probably less interesting. In theories with fc-dependent 
linear growth, only the peaks with exactly the right large 
scale surroundings (determined by the fc-dependence of 
linear theory) will survive at later times; this raises the 
possibility that the correlation between galaxy clusters 
and their environments can constrain theories of large- 
scale modifications to gravity. We have not pursued this 
further, but note that this is consistent with recent anal- 
yses of dark matter halos 



To conclude, it is worth mentioning that halo-based 
approaches, which provide a reasonably good description 
of the weakly nonlinear clustering of simulated dark mat- 
ter haloes and galaxies [e.g. HI], are now commonly used 
to extract cosmological information from redshift surveys 
(see [H| for a review). Although the dark matter halos 
are the local density maxima of the evolved matter dis- 
tribution, there is no easy correspondence with the ini- 
tial density maxima. This is the reason why the peak 
model has somewhat fallen out of favour. We believe our 
work has shown that many insi ghts can be gained from a 
study of density peaks [e.g. [13, 113], particularly with re- 
gard to a number of effects - including scale dependence 
and stochasticity of the spatial and velocity bias - which 
matter in the age of precision cosmology. 
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Appendix A: Checking the consistency of the peak 
biasing relation 

In this Appendix, we sketch the derivation of the cross- 
correlation between peaks and the underlying density 
field, £ p k,s{r), and the averaged peak pairwise velocity, 
V12 (r, /i) , which have not been derived previously. We 
compute both quantities using the peak constraint and 
demonstrate that the results are consistent with those in- 
ferred (after a trivial calculation) from the peak biasing 
relation See [29| for complementary details about 

the calculation. 



Cross-correlation between peaks and density 
field 



Let t]i — did(x)/cr 1 and = didjS(x)/(T2 be the 
normalised first and second derivatives of the density 
field. Furthermore, let A be the diagonal matrix of en- 
tries diag(Ai, A2, A3) where Ai > A2 > A3 is the non- 
increasing sequence of eigenvalues of the symmetric ma- 
trix — £. The cross-correlation ^ p k, 5 (r) = (5n p k(xi)v(xn)) 
(r = \x2 — Xi\) follows from the Kac-Rice formula |90l |. 
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where n p k is the differential number density of peaks of 
height v as given in Eq. ©, R\ is the typical radius of 
peaks and A3 is the smallest eigenvalue of — C a t x\. We 
have omitted the ^-dependence for brevity. The Heav- 
iside step-function arises because we are interested in 
counting maxima solely, for which Cii is negative defi- 
nite at the extrema position. Unlike [2(| who expresses 
the covariance matrix in a coordinate system where the 
two density maxima lie on the z-axis, we write the expec- 
tation value in the right-hand side of Eq. (IA1[) as an inte- 
gral over the angular average, joint probability distribu- 
tion function P(y 1 v 2 ;r). Here, y T — Ca) is a ten- 
dimensional vector whose components £a> A — 1, . . . , 6 
symbolise the entries ij = 11, 22, 33, 12, 13, 23 of Qj. 

To evaluate the 11-dimensional covariance matrix 
C(r), it is convenient to split the degrees of freedom asso- 
ciated with the tensor ( into the scalar u = — tr£ = J2 i \ 
(so that u is positive when A3 > 0) and the traceless ma- 
trix £ = £ — l/3(tr£)I , where I is the 3x3 identity 



matrix. Let (a designate the 5 degrees of freedom of £. 
We see that C(r) = (l/47r) JdQ- r C(r) has the block diag- 
onal decomposition C = diag(Ci, C2, C3). Consequently, 
P(Vi, v 2l r) can be expressed as a product of three joint 
probability distributions 



P(y, v 2 \r) = P{v u uy,V2\ r) P( m ) P(Ci) 



(A2) 



where, for shorthand convenience, subscripts denote the 
variables evaluated at different Lagrangian positions. 
The 1-point probability distributions P(rji) and P(Ci) 
are 



P(m) = 



3/2 



exp 



2 



(A3) 
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while the joint density P{vi,Ui,V2]r) has a covariance 
matrix 



Ci = 



1 



(i). 



1 

7i 



7i 
1 



(A4) 



Upon inversion of Ci , the quadratic form Qi that appears 
in the probability density P{v\, u\, i>2] r) reads as 



where 



2(1-7?) 



2A* 



At = 1- 



1 



AO) 
so 



7i So 



Ax = — 



1 



' 1 

.°o 

1 

(T 2 



l-ll 

v\ - 7i«i 



1 



Ml 



"7i 
71^1 



So 



(0) 
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(A6) 
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The calculation now proceeds along lines similar to [29] ]. 

We choose a coordinate system whose axes are aligned 
with the principal frame of £1 and introduce the asym- 
metry parameters 



« = (Ai-Aa)/2 

w = (A x - 2A 2 + A 3 ) /2 



(A8) 



Our choice of ordering impose the constraints v > and 
— z < w < v, while the peak constraint enforces (u + 
w) > 3w. Upon integration over the angular variables 
that define the orientation of the orthonormal triad of £1 
and the variables v and w, the cross-correlation £ p k,<5(f) 
of peaks of height v is given by 



= cr G (7l,7l^l) 
p 00 

duif(ui) 



+00 



dv 2 v 2 - 



-(u 2 -A x ) 2 /2A e 



e -(u 1 -7l^l) 2 /2(l-7 1 2 ) 

^(1-7?) 



(A9) 
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where the auxiliary function f{u) is defined as in 
Eq. (A15) of [H, and 

r°° e -(z-^) 2 /2(i-7 2 ) 
G„( 7 ,w)= / dxx n f(x) — —=—-. (A10) 
Jo A/27r(l-7 2 ) 

are moments of the peak curvature. In particular, u{v) = 
Gi (71 , 71 v) /Go (71 , 71 1 / ) is the average curvature of peaks 
of height v. Finally, the integral over v 2 is performed and 
we arrive at the desired result: 



Cpk.aO) 



1 {y - 7ii2) (0) 1 (jj- 7lI/ ) (1) 

a (1 - 7 f) *» W (72 (1 7?) f ° ( J • 

(All) 

This agrees with Eq. (IT51) , which was obtained with much 
less effort from the peak biasing relation (|4]). It is worth 
noticing that, while the derivation based on the peak 
biasing relation is formally exact at the first order only, 
this appendix shows that Eq. (jAllI) is exact to all orders. 

The cross-correlation £ p k,<5(v) has a straightforward in- 
terpretation: it is the average density profile around den- 
sity maxima, i.e. (5(a;2)|peak at cci). As shown by [l7j . 
this constrained density profile can be calculated eas- 
ily and, after some algebra, one indeed finds Cpk.aM = 
(£(a: 2 )|peak at x x ). Note that ip(r) in Eq. (7.10) of ref- 
erence [I?} corresponds to our Cq ^?"). Therefore, their 
Eq. (7.10) appears to have an additional factor of 1/3 

which multiplies the factors of V 2 ^ ' = — £g in our 
expression - but this is only because they measure r in 
units of R\ - there is, in fact, no difference. 



2. Mean streaming of peak pairs 

The calculation of the mean streaming is more involved 
since we have three more degrees of freedom and an extra 
angular dependence. Our derivation is based on refer- 
ence , who calculated the pairwise velocity dispersion 
along the line of sight. 

We introduce the normalised velocity field wi = 
Nij (aH f<J-i) = Vij(j-\. Also, we assume that the line of 
sight axis coincides with the third-axis, such that Aw z 
denotes the difference 1373(332) — voj,{x\). The line of sight 



pairwise velocity weighted over all peak pairs with co- 
moving separation r can be expressed as 



[1 + Cpk(f)] v 12 (r, n) = n pk 2 cr_! 



(A12) 
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>dy 1 dy 2 Am z n p]i (x 1 )np k (x 2 )P(y 1 , y 2 ; r) 



where /1 is the cosine of the angle between r = r/r 
and the third axis, and <f> is the azimuthal angle in the 
plane perpendicular to the line of sight. The local peak 
density n pk (x) is given by 3 3 ^ 2 Ri 3 \det((x)\S 3 [ri(x)}, sup- 
plemented by the appropriate conditions to select those 
maxima with a certain threshold height. 

In the above expression, the joint probability density 

-P(yi)2/2; r ) is now a function of y T = (u7i,rft,i>,(A), 
where rji = owing to the peak constraint. The cor- 
responding covariance matrix C can be decomposed into 
four 13 x 13 block matrices, the zero-point contribution 
M in the top left and bottom right corners, and the cross- 
correlation matrix B(r) and its transpose in the bottom 
left and top right corners, respectively. In the large dis- 
tance limit r» 1 where |B| <C M, an expansion in the 
small perturbation M yields at first order 



(2 7 r) 13 |detC| 1 /2 



(l + y7M- 1 BM- 1 y 2 ) 
(A13) 



where the quadratic form Q(y 1 , y 2 ) reads 
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3^ 2 



(71Z/1 + trCi)' 
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vj\ being the velocity vector at comoving position x\. 
The inverse M _1 and B(r) can be further decomposed 
into the block matrices 



M = 



P R [ 
R Q 



B = ( Bl 

B3 B2 



(A15) 



where 
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The explicit expressions for are too long to be given here as they depend upon the correlation functions of Wi, 
rji, v and (,A in a rather complicated way. Fortunately, the mean streaming involves only the azimuthal average 
B(r, (j.) = 1/(2-71") Jd(f>B(r), which can generally be expanded as 



B(r )M )=^^(r)L,(M), &(r) = 



B 



B| T 



T>£ Tit 



(A18) 



Note that the multipole matrices B| and B| T are not independent of each other since we have B| = (— l^B^. As 
we will see shortly, all the contributions but that from the £ — 1 multipole (unsurprisingly) cancel out. We will thus 
detail the results for the dipole contribution solely. After some algebra, we find 
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whereas the matrix B3, is identically zero. Right and left multiplication by M 1 then gives 
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where the functions oti{r) are identical to those defined in Eq. (A12) of [29[. Namely, 
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The rest of the calculation is easily accomplished owing 
to the separability of the one-point probability distribu- 
tion P(y) into the product P^(vJi)P u c,(v, (a), where P v q 
is the one-point distribution of the density and its sec- 
ond derivatives (The first derivatives merely contributes 
a normalisation factor), and 



33/2 



(2 7 r)3/2(l- 7 2 ) 3 / 2 



exp 



2(l-7o 2 ) 



(A22) 

is the velocity distribution of density peaks. In particu- 
lar, the first moment vanishes while the second moment 
(zuf) is the one-dimensional velocity dispersion of density 
maxima, 



K 2 ) 
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7o 2 ) 



(A23) 



The scalar yjM~ 1 BM~ 1 y 2 contains terms linear and 
quadratic in Wi as well as terms independent of the ve- 
locity. Upon multiplication by Arr7 z and integration over 
the velocities, only quadratic terms survive. We eventu- 
ally find 

Jd^vj^^Amz (yjM^BM^y^j P ro (^i)P ro (^ 2 ) 

= [ai (vi + t> 2 ) + 7i"i ( tr Ci + t< 2 ) - 3a 3 (£1,3 + Cz,3)] 
x (l- 7o 2 ) Li(/x) 
3 



[Ci,i + C 2 ,i + Ci, 2 + C 2 , 2 - 2 (Cl,3 + C2, 3 )] 



.(1/2) _ (7= (3/2) 
?3 S??3 



(A24) 



CT-102 

Here, Ci.A arL d C, 2 ,a designate the component (^a of the 
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hessian £ at location Xi and respectively. As we can 
see, although the even multipoles cancel out, a term pro- 
portional to Lz{n) remains. Also, the dipole receives a 
contribution from —303(^1^ + ^2,3) which is not invariant 
under rotations. However, we have to remember that the 
principal axes of the tensors £1 = Ci x i) an d C2 = C( a; 2) 
are not necessarily aligned with those of the coordinate 
frame. Let us first consider £1. Without loss of general- 
ity, we can write £1 = — OAO T , where O is an orthogo- 
nal matrix and A is the diagonal matrix consisting of the 
three ordered eigenvalues \ of — (1 . The properties of the 
trace implies that tr£i = — trA, while dj = — Yli 
Since the one-point probability density P(y) does not 
depend on O, the integral over the SO(3) manifold that 
describes the orientation of the orthonormal triad of £1 
is immediate, 



3 3 
J i=l J 6 i=l 



1 



rtrCi 



(A25) 

Similarly, averaging over the orientation of the eigenvec- 
tors of C2 yields fdO&j = (l/3)tr£ 2 - Consequently, the 
£ = 3 term vanishes and we only need to integrate 



ot\ {v\ + v 2 ) + (701 - a 3 ) (trCi + trC 2 ) (l - To) ^i(m) 

(A26) 

over the eigenvalues of fi and £2 subjects to the peak 
constraint. Substituting the expressions ([5]) of the bias 
parameters b v and b^, the result can be recast into the 
form of Eq. (1341) when v\ = v<i = v. 



Appendix B: The cross-correlation coefficient for 
powerlaw spectra 

In this Appendix, we examine the large-scale behaviour 
of the cross-correlation coefficient for density peaks as- 



suming a power law spectrum of density fluctuations. 

The cross-correlation coefficient in configuration 
space can be written as 



|(r) 



1 + H(r) ' 



TZ(r) = 



t(i) /tf (o) 



bu/h 



(Bl) 



It is larger than unity when TZ < 0, i.e. when £^ /(.^ < - 
[Co^/Cq ^ 2 - In Cold Dark Matter cosmolo gies, the cor- 
relation functions £q must be calculated numerically 
because the spectral index is a smooth function of 
wavenumber. For a no-wiggle powerlaw power spectrum 
P$(k) = A s k n3 however, they take the exact form 



ft 



(B2) 



where T(a) and ii*i(a, 7; — 2) are the Gamma and con- 
fluent hypergeomctric function in the arguments a — 
n + 3/2 + n s /2, 7 = 3/2, and z = r 2 /(4i?|); and i? s 
is the characteristic radius of the window function as- 
sumed Gaussian. Since Eq. (IB1|) only holds at large 
separation r» 1, we consider the limit \z\ — > 00 to the 
above expression, in which i-Fi(a, 7; —z) has the follow- 
ing asymptotic expansion (in a suitable domain of the 
complex plane |91j). 



IY7) 
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k=l 



T(fc + a)r(fc + a-7 + l) _ fc 
fc!r(a)T(a-7 + 1) Z 



In the right half-plane of the variable z (i.e. for Re(z) > 
0), the first term in the right-hand side of Eq. (|B3|) is the 
subordinate part of the asymptotics, whereas the second 
is the dominant part. We will now calculate 1Z(r) for 
a few integer values of the spectral index covering the 
range [—3, 2] shown in Fig. [7] As we will see shortly, the 
cross-correlation coefficient can exceed unity even if the 
underlying power spectrum P(k) is a featureless power 



law. The exact amount of stochasticity, however, crit- 
ically depends upon the shape of the underlying power 
spectrum. 

In the particular case of a white noise spectrum, n s = 

0, the dominant part cancels out owing to the fact that 
r(7 — a) has simple poles at 7 — a = n = 0, 1, 2, • • • , 

1. e. r(-n)- 1 = 0. Moreover, T(k + 7 - a)/T(-f - a) = 
(k — 1 + 7 — a)(k — 2 + 7 — a) x ■ ■ • x (7 — a) vanishes 
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when k > n + 1, so the summation in the subordinate 
part involves a few terms solely. In fact, the asymptotic 
expansion gives the exact result, 



87T 3 / 2 i?| 
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which yields 



K(r) 



b»R% , 3 
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As can be seen, lZ(r) becomes negative at separation r > 
VSRs, so the cross-correlation coefficient is greater than 
unity at large scales. Note, however, that the dominant 
part is nonzero for any small n s different from zero. More 
precisely, upon writing n s — e where < |ej <C 1 and 
momentarily ignoring a factor of A s j (47r 2 )i?^ 2a , we have 



r (n 



r (-»-§) 

The sign of TZ(r) is equal to that of 



r (3/ 2 ) ,-n-i 
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On inserting these expressions into Eq. (|BI[) . we find 
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Again, r% > I at sufficiently large separation r 3> 1. 
Note, however, that 1Z decays much more rapidly to zero 
when n a = —2. Furthermore, one can show that < 1 
for n s = —2 + e, and > 1 for n s — —2 — e, where 
0<e<l. In other words, there is a jump discontinuity 
at n s = —2. 

When the spectral index is an odd integer, e.g. n s = 
— 3,±1, the subordinate, complex-valued part is expo- 
nentially suppressed relative to the dominant, real- valued 
part. For n s = — 1, we find 
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upon including the first two terms of the dominant part. 
After some simplification, we arrive at 



(B13) 



which is positive for any small nonzero e. Therefore, the 
cross-correlation coefficient is discontinuous at n s = 0. 
The same analysis also shows there is a similar disconti- 
nuity at n s — 2. These discontinuity points are marked 
as empty symbols in Fig. [7] 

For n s = —2, the dominant part is non- vanishing only 
when n = 0. Furthermore, for n = 1 and 2, the subordi- 
nate part only sums a finite number of terms. Explicitly, 
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Similarly, we obtain 
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(B14) 



for n., = —1. In both cases, 7?. > so the cross-correlation 
coefficient is less than unity at large scales. Finally, for 
n a = —3, the density correlation Q (r) diverges owing 
to the presence of T(a) = r(n). Consequently, the cross- 
correlation coefficient is unity at all scales. 
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